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In this paper we derive a general expression for the transmission coefficient using the method 
of reactive flux for a particle coupled to a harmonic bath surmounting a one dimensional inverted 
parabolic barrier. Unlike Kohen and Tannor [J. Chem. Phys. 103, 6013 (1995)] we use a normal 
mode analysis where the unstable and the other modes have a complete physical meaning. Impor- 
tantly our approach results a very general expression for the time dependent transmission coefficient 
not restricted to overdamped limit. Once the spectral density for the problem is known one can use 
our formula to evaluate the time dependent transmission coefficient. We have done the calculations 
with time dependent friction used by Xie [Phys. Rev. Lett 93, 180603 (2004)] and also the one used 
by Kohen and Tannor [J. Chem. Phys. 103, 6013 (1995)]. Like the formula of Kohen and Tannor 
our formula also reproduces the results of transition state theory as well as the Kramers theory in 
the limits t — > and t — > oo respectively. 



I. INTRODUCTION 



The dynamics of a particle surmounting a barrier is 
an important and interesting problem in chemistry and 
physics. In chemistry, chemical reactions are most com- 
mon examples where one studies the dynamics of bar- 
rier crossing along a suitably defined reaction coordinate. 
While in physics nucleation phenomena, electrical trans- 
port etc. involve barrier crossing. A review by Hanggi, 
Talkner and Borkovec [l[ gives a comprehensive account 
on the subject of barrier crossing. The simplest possible 
approach to calculate the rate constant for such processes 
is given by transition state theory [H, . As the simplest 
version of the problem one can describe the rate pro- 
cess along a single reaction co-ordinate. Then the one 
dimensional version of the transition state rate constant 

krsT reads as, kxsr = f^ ex P[ — fcfrb wnere w o is the 
angular frequency in the reactant well and is the bar- 
rier height or the activation energy. The derivation of 
transition state rate constant is based on two assump- 
tions. Firstly an overall thermal equilibrium is assumed 
and secondly frictional effect is not taken into account by 
incorporating what is known as "no recrossing" assump- 
tion. Subsequently, Kramers [l|, 0, [H, 0] theory takes into 
account of the friction in a phenomenological way. His 
formulation was based on Markovian dynamics of a Brow- 
nian particle escaping from a metastable state. Further 
extensions of Kramers's theory were carried out by Grote 
and Hynes [5| followed by Hanggi and Mojtabi [6[. They 
took into account of the non-Markovian nature of the dy- 
namics. Another approach to calculate the rate constant 
is due to Chandler [3, [1] known as the method of reactive 
flux in chemical physics literature. In this approach the 
thermal rate constant is written as a correlation function 
which is nothing but the ensemble average over infinite 
number of trajectories starting at the barrier top and 
ending on the product side at time t. The rate constant 
at time t is written as a product of the time dependent 
transmission coefficient n(t) and the transition state rate 



constant krsr- Thus one writes, k(t) = n(t)kTST- hi the 
limit t — > 0, k{£) — > 1 and one obtains fc(0) — kxsr- Here 
we follow the method of reactive flux to calculate the time 
dependent rate constant for a particle bilinearly coupled 
to a harmonic bath Q [rJ [n|. The bilinear coupling 
allows us to get an analytical expression for the time de- 
pendent transmission coefficient. Although the method 
of reactive flux and the hamiltonian we have used are well 
known in the literature, as far as our knowledge nobody 
has used such hamiltonian to calculate the time depen- 
dent transmission coefficient analytically. The calcula- 
tion is of interest as the transmission coefficient for the 
overdamped limit of the model was obtained recently by 
simulation [12]. Our approach directly gives expression 
for the time dependent transmission coefficient in terms 
of the spectral density J(ui) and not restricted to over- 
damped limit. A normal mode analysis of the coupled 
Hamiltonian was made by Pollak (l3l . [Til HH and others 
to analyze the Kramers problem. Here we point out that 
the approximation is very powerful and can be used to 
get the exact expression for n(t). Further we believe that 
it can also be used to get the quantum expression for 
n(t). An exact calculation of K(t) using the phase space 
distribution function formulation was carried out in clas- 
sical regime by Kohen and Tannor [lj| |T3] , Bao , and 
more recently in the context of single enzyme kinetics by 
Chowdhury and Cherayil [l9( . Also the calculation of the 
transmission coefficient in the quantum domain is carried 
out by Ray's group [201 ] . Our analysis in terms of normal 
modes has the advantage that it makes the physical ideas 
very clear. The paper is arranged as follows. In section 
II we introduce our system plus harmonic bath hamilto- 
nian and the Generalized Langevin Equation (GLE). In 
section III we mention the method of reactive flux and 
using our hamiltonian derive an analytical expression for 
the time dependent transmission coefficient, n(t). Sec- 
tion IV deals with the calculation of n(t) for Markovian 
and Non Markovian dynamics. Section V is conclusion. 
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II. GENERALIZED LANGEVIN EQUATION 
AND NORMAL MODE ANALYSIS 

We consider a particle coupled to a harmonic bath 0, 
Hoj with the total hamiltonian 



P 2 \ N ( P 2 

3=1 V J 



+ mjuj? ( 



— <3 



rrijUj 



Here (Pj, Qj) are the momenta and coordinates of the 
jth bath oscillator whose mass and frequency are rrij , ir- 
respectively. (P, Q) is the momentum and coordinate of 
the system, cj couples the bath oscillator to the system. 
Now in terms of the mass weighted coordinates 



q = VMQ, qj = ^/mjQj 
our Hamiltonian becomes 



(2) 



P 



N 




and then with this Hamiltonian one can derive what 
is know as the Generalized Langevin Equation (GLE) 
which reads as 



Kt) + V\q) + I dt' 7 (t - t)q(t) = C(t) (4) 



One can also write the time dependent friction in terms 
of the spectral density of the bath J(ui) as follows 
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where 



where U_ is an orthogonal matrix such that U_D.IL is 
diagonal. D is a dynamical matrix given by 
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The reaction co-ordinate q may then be written as 

3=0 

Of the modes rjj, we take r?o to be the unstable mode. 
We denote the frequency of rjj as Xj . rj has the imaginary 
frequency and we write it as Ao = iA, where A is real. As 
rjj is a normal mode, its time development is given by 



Vj(t) =Vo(0)cos(X j t)+p rij (0) 



sin(Ajf) 
A., 



with j = 0,1,2, N. 

Now the Hamiltonian when written in terms of the 
normal modes the coupling is no more there and is given 
by 




(9) 



3=0 



In the next section we analytically derive the rate con- 
stant in the normal mode description using the above 
Hamiltonian. 



III. METHOD OF REACTIVE FLUX AND THE 
RATE CONSTANT 

In this section we briefly describe the method of reac- 
tive flux @,Q to calculate the rate constant. The method 
of reactive flux expresses the thermal rate constant (writ- 
ten in the mass-weighted coordinates) as 



N 



2 ^ mjWj 



5(uj - uj) 



One can easily show that 



(C(*i)C(*a)) = MksT^h 



(6) 



(7) 



We take the barrier top to be at q — and introduce 
lu 2 = (^jr^j ■ For small amplitude motion around the 

N 

barrier top, one can use the normal modes, T)j = Ujiqi 

i=o 



k(t) 



(q(t)S(q(O))0p(q(t))) _ k n (t) 

k d (t) 



(0«(<?(o))> 



(10) 



where the angular bracket indicates the thermal av- 
erage over the system as well as the bath degrees of 
freedom. The top of the barrier is chosen as q(Q) = 0. 
9n(q(0)) is 1 if q(0) < 1 and otherwise. 9 p is just 1-9 r, 
q is the velocity of the particle. 

The numerator of Eq. (TlT)|) could be written as 



/.■„(/) = ( ^(exp(ir q(0))6(q(0))6 P (q(t))) 



r =0 
(11) 
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Expressing the step function as an integral over delta delta function the expression for k„(t) becomes 
function and using the Fourier integral representation of 

I 



oo oo 



k n (t) = 



(27T) 



— 2 / dh(—i)-^— I dso I dsi exp(— is±h) (exp(ir q(0) + isiq(t) + is o q(0))) 



dr a 



— oo — oo 



ro=0 
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The angular bracket indicates the thermal average over exp(—/3H nor ) , and carrying out the integration over all 
all degrees of freedom of the problem and it is done by coordinates and momenta degrees of freedom. Writing 
multiplying the quantities within the angular bracket by k n (t) in terms of rjj leads to 



oo oo 



kn(t) 
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where is the barrier height. After one carries out the integrations the expression for k n (t) can be written as 
I 



JV 
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Similarly the denominator of the Eq. (fit)]) can be writ- ten as 
I 



*d(*) = n y ^y ^ ex p(^(y + i-^)) = Il(^f 

i — 0— no — oo ^ ' J — \ / 



(14) 



Then the time dependent rate constant becomes 



kit) 



N 

k n (t) = 1 /=Q 3 

fc d (f) 2ttA 



-c(f) 



g Aj V(c(*) 2 -c(o)^ 

i=i 



T exp(-/3 J B b ) 



(15) 



4 



n Co A? 

Here we use the following identity [l|, — = 

a n h 

3 = 1 

and this simplifies our k(t) 



In the next section we calculate c(t) and from that n(t) 
for Markovian as well as for Non Markovian noise. 



-m 



2lTUJ b V (c(i) 2 - C(0) 2 ) 



r eM~(3E b ) (16) 



where, c(t) = I £ ^ c °f* f) 



As usual [16[ is 



v 3=0 



defined by 



fc(t) = K(t)fc(0) 



(17) 



We show in the Appendix that fc(0) = ^ exp(— /315&), 
which is nothing but the transition state rate constant 
krST- In the Appendix we also show how one gets 
Kramers rate constant in the t — > oo limit. 

Using the above expression for fc(0) our time dependent 
transmission coefficient, n(t) becomes 



1 



«(*) = — 



-c(t) 



V(c(t)2- C (0)2)2 



(18) 
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IV. CALCULATION OF c(t) AND /s(t) 

Calculation of nit) involves the calculation of c(t). 
First we describe a general derivation of c(t) then we 
evaluate nit) for different cases. 



A. Calculation of c(t) 



c{t) is defined by c{t) — ^2 Uj ° °°^ Ajt ^ . One proceeds 



N 



N 



by writing £ [/ 2 <5(u, 2 - A 2 ) = -ilm £ 



To 



JV 



here we adopt notations similar 



evaluate £ ? 2 , J ' , 2 \ 

to quantum mechanics and also use the fact that D is a 
matrix with eigen values A 2 , j = 1,2, AT. 



JV 



- (w* + i V -tf ^ (0|J> (c 2 + z ?/ -A 2 ) 



(i|o) = ^<o 
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(oj 2 + it] — D) 



(J I 0) 



Then by introducing |j) (j| in between the sum and 



using the resolution of identity ^ \j) (j\ = I. 

3=1 



(19) 



j:^ 2 -a 3 2 ) 

3=1 
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(lu 2 + ir/ — D) 
I 



|0> = — Im (Q(uj 2 + irfi) 



(20) 



get 



Using partitioning technique one can evaluate (g_) Q0 to Ag __i Im (G(cj 2 )) = |] U\5{uj 2 - A 2 ) = p(w 2 ) We 



(10 



- 2 + - b 2 -i^E 



(21) 



^ mjw|(w 2 +ir;-w?) 



Then using the definition of J(w) Eq.© it becomes 
& 2 +«?))„o 



+ UJ 



2 _ 2uP_ r 

b Mir J uj(u 2 +iri-Cj 2 ) 




(22) 



3=0 

OO 

can write c(f) = 2 / ^ cos(^)p(w 2 ). 
o 



B. Calculation of n(t) with T (t) = y2H (2H - 1) |i| 2iM 

First we perform the calculation with the time depen- 
dent friction j(t) used by Xie [2l[, which is 
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FIG. 1: Time dependent transmission coefficient, «(t) against 
time, for H = \ (dashed line) and H = | (solid line), param- 
eters used are Lib = 1,7 = 0.1 



a inverse Fourier transform of Eq. ([5]) and write J(u>) in 
terms of 7(i). 



(24) 



J(u) = Moj J dtj (t) cos(ut) 
o 

Using the definition of ^(t) we get 



J(uj)duj 



2uj 2 



Mir J uj (uj 2 + if] — uj 2 ) 
o 



-yT{2H + 1)A 



2-2H 



First we consider the case a-> 2 = —A 2 which means we 
consider the unstable mode of the problem. Then 




t 



FIG. 2: Time dependent transmission coefficient, K,(i) against 
time for Non-Markovian friction, 7(t) = yhexp ( — ht). Three 
different regimes, caging regime (solid line) with parameters, 
7 = 150, h = 0.01, uib = 1, intermediate regime (dotted line) 
with parameters 7 = 5, h = l,u>b = 1, and nonadiabatic 
regime (dashed line) with parameters 7 = 50, h = 0.01, u)b = 
1. 

7 (i) =j2H(2H -l)\t\ 2H ~ 2 (23) 

with, i < H < 1, thus the thermal fluctuations are de- 
scribed by fractional Gaussian noise [2l[. Now we take 



(0- a2 ))oo = (- a2 + "2 - ^ 2H + i^-^r 1 

(25) 

Our next task is to find out the poles, i.e. the unstable 
modes for which G 00 (—A 2 ) blows up, i.e. the solutions 
of the equation 

- A 2 + uj 2 - jT(2H + 1)A 2 - 2H = (26) 

On the other hand, to find p(uj 2 ) we have to consider a 
situation where ui 2 is positive. In other words the stable 

modes. Now the Evaluation of p(uj 2 ) involves the inte- 

00 

S ral ]h£ / Usin S the identity ( „ 2+ >_ aa) = 

P ( jjihfii ) -inS (v 2 _ <^ 2 ) > where P f ^^rzg? ) is the prin- 
cipal value of ( ) we find 



00 

2 4- I -2, = - «P r(2H+2). 2 - 2H 1 

Mir J uj {uj 2 + irj - uj 2 ) y ' v ; 



(27) 

Thus 

& 2 ))oo = + LU b+ CX P {iHir) T(2H + 2)uj 2 - 2H 1 y 1 

(28) 

and 



p(uj 2 ) = -±Im(g(uj 2 )) 



on 



7r(2g+l)c^ J 



-l(Hit) 



(29) 



^ ((w 2 +w2 +7 r(2H+l)w 2 - 2 H cos(ff7r)) 2 +7 2 sin 2 (Hn)r(2H+l) 2 uj 4 - 4H J 
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1. The Case H=\ 

First we consider H — |, for which the j(t) = 2"f6(t) 
and with this the noise-noise correlation function defined 
in Eq. ([7]) becomes delta function correlated, reading 



(C(*i)C(*a)) = 2Mjk B TS(h - h) 



(30) 



Then the dynamics is described by usual Langevin equa- 
tion rather than a GLE. Thus the dynamics is then 
Markovian i.e. with no memory. For the Markovian 



case Eq. (1261) becomes, — A A 



7 A = and has 



two roots. But only one of them is physically accept- 
able. The root is, A, 



^ = -i + V(i) 

from Eq. ([29]) one gets p(uj 2 ) H= i 



uji . Similarly 



Then A ff= .i an d p( uj2 )h=± are used to calculate c(t) and 
taking t — > limit one gets c(0). In this case we could 
do these analytically. The time dependent transmission 
coefficient in this case is given by 



K(t) 



(_l + e t</-7)) 



yj(l - 4ef + 2e*(/-T) + 2e 2t (f-^) - (^)(/ - fe 2 ^-^ + 2 7 (- 



s /t)) 



1/2 

where / = 7 + {Auj b 2 + 7) . 

n(t) is then calculated using the above expression and 
plotted against time (Fig. I). 



2. The Case H 



When H ^ | then the noise is no longer delta function 
correlated, i.e. not a white noise with Gaussian distri- 
bution. This regime is know as Non Markovian regime 
which has memory. We have similarly evaluated K,(t) for 
H = I and plotted against time in Fig. 2. But in this 
case the integrations are done numerically since analyti- 
cal evaluation of integrals were not possible. 



C. Calculation of ft(t) with 7 (£) = *yhexp (—ht) 



In this section we carry out the calculation of n(t) with 
a time dependent friction j(t) — jhexy (—ht). The same 
friction is used by Kohen and Tannor [16|, [l7| . First we 
calculate J(uj) with this friction as follows 



J(uj) = M7/1W / exp (—ht) cos (ujt) dt = 



■yMh 2 uj 
(h 2 +u> 2 ) 



Then the following integral is evaluated as follows 



2uj f J(uj)duj 2uj h f duo ^huj (uj — ih) 



Mir J (a 2 + uj 2 ) (uj 2 + ir\ - uj 2 ) tt J (a 2 + uj 2 ) (uj 2 + ir\ - uj 2 ) (h 2 + uj 2 ) 




We have used the identity, = P ( ) - and 

5 (uj 2 — Co 2 ) , where P ( rjTLZjp ) is the principal value of 



ITTi 

Thus 



(w 2 -lD 2 )- 



t nl 2x\ (2,2 ihio (cj - ih) \ 1 



, I , , „ NN Lo-ih 2 (h 2 + uj 2 ) 

p(uj 2 ) = --Im(g(0J 2 )) c - 



'"" TT ({{h 2 + UJ 2 ) (h 2 + UJ 2 ) - -fhuj 2 ) 2 + UJ 2 -f 2 h^ 



One should notice that 
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lim 



UJ-fh 2 (h 2 + LO 2 ) 



TT (((/l 2 + LO 2 ) (h 2 + LO 2 ) - -fhu 2 ) 2 + L0 2 J 2 h^ T ((W 2 + L0 2 b ) + J 2 UJ 2 ) 



Thus in this limit the result obtained is identical with 
the friction y(t) = j2H (2H - 1) \t\ 2H ~ 2 with H = ±. In 
other words the Markovian dynamics is recovered. 

To calculate n(t) we proceed as in the previous case. 

First considering the case lo 2 = —A 2 , one gets 

2uj 2 J J{uj)dw _ 7/1A 

~Mtt J (a 2 + Cj 2 ) {uj 2 + it] - Co 2 ) ~ (h + A) 
o 

Thus 

Hence to find out the poles one has to solve the solution 
of the equation 



method of reactive flux 0, Q to calculate the time de- 
pendent transmission coefficient, n(t). Analytically we 
derive a general formula for n(t) Eq. (TTg]) . For the 
Markovian case with no memory (H = i) for the fric- 
tion j(t) = j2H(2H-l) \t\ 2H ~ 2 with H = \ which is 
same as with h = oo for the friction j(t) = 7/1 exp (—hi), 
n{t) is calculated analytically. Whereas in the case with 
memory {H = |) and with any finite value of h the cal- 
culations were done numerically. In all the cases a plot 
of n{t) vs t starts from 1 at t = and reaches a plateau 
in the long time limit. As expected our formula for the 
time dependent rate constant, k(t) becomes equal to the 
transition state rate constant when one takes the limit 
( -> 0, 0~7| ■ Similarly Kramers rate constant is ob- 
tained by taking t — ► 00 [16L [TtJ • In future we would like 
to extend our formulation to quantum domain. 



- A 



7/1A 
(h + A) 







(31) 



To explore the non-Markovian dynamics we choose 
three set of parameters corresponding to three different 
regimes of non-Markovian dynamics following Kohen and 
Tannor . fl7|. These regimes are called non-adiabatic, 
caging and intermediate according to Kohen and Tan- 
nor [3, The time dependent coefficients in these 
regimes are calculated numerically and plotted against 
time. We observe similar kind of oscillatory behavior of 
n(t) in the caging regime as obtained by them. Similarly 
in the intermediate regime the plot looks similar with a 
nonmonotonic decay (Fig. 2). 



CONCLUSIONS 
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VII. APPENDIX : TRANSITION STATE RATE 
AND KRAMERS RATE FROM k(t) 



The paper shows an elegant way of combining the 
traditional system plus reservoir model 0, [l(3| and the 



First we take the limit t — > 0. Then from Eq. (jT6|) one 
writes 
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*(0) = Hm^o fc(Q + 1) = lim^o (e(0+ ^ )a)A «p(- W = 
Hm^o ^ exp(-m) 

^(c(0)+tc(0) + ^c(0)) -c(0)=J 

i^limt^o zmi r eM-PE b ) = j-^a exp(-8E b ) = — exp(~-(3E b ) 



Where we have used the facts, c(O) = and c(O) = — 1 i,,2 _ , ¥ , , m\ — i 

— u jQ — — i. inow one can also snow, — c^uj — ~ 



This is because c(t) = —2 J dujp{uj 2 ) sin(a;i) making 



o 



4=0 
as follows 



c(0) = 0. Similarly c(t) = — 2 J dtu p(tu 2 ) sm(cut)ui 

o 



N 

c(0) = ^oA- 2 c/ 0j = E (o I j) A7 2 (J I o) 

4=0 J 3=0 4=0 

W JV 

lim^^o E (0 I 3) 01 T^uy |i) (j I 0) = - lim^ £ (0| |0) = - lim w ^o (£(w 2 )) 00 

4=0 j=0 



I 

From Eq. (|22|) one can see that lim w 2^ (G(w 2 )) 00 = Similarly taking t — > 00 one gets the Kramers rate 

4r = -c(0). Then fc(0) is nothing but the transition constant as follows 
state rate constant, Utst- 

fc(0) = ^ exp(-/3£ b ) = k TST (32) 



fc(oo) 



1 <^0 
27T 0->(, 



ah(At) 



^^ cosh( A t) y 



exp(-/3-E 6 



£^(Atanh(Ai)) t 



exp(-/3£ , i) ) 



Tnen Now if one replaces A by A H= i = —2- + \f{^)' 

J_^0 

27T Wfo 



fc(oo) = A—— ex P (-« (i/ = 3 meanS White n0iS6) ° ne g6tS 



fe(oo) = -? + V (i) + £^«p(-/32&) 



Which is nothing but the usual Kramers rate. Also in 
the overdamped limit i.e. when 7 >> ui 2 one gets 
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«„, _ (_2 + 2 ( 1 + i t (£))) ^e*p<-« 

I 



which is Kramers rate expression in overdamped 
regime. 
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